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ABSTRACT 

Context. In some galaxy clusters powerful AGN have blown bubbles with cluster scale extent into the ambient medium. The main pressure support 
of these bubbles is not known to date, but cosmic rays are a viable possibility. For such a scenario copious gamma-ray emission is expected as a 
tracer of cosmic rays from these systems. 

Aims. Hydra A, the closest galaxy cluster hosting a cluster scale AGN outburst, located at a redshift of 0.0538, is investigated for being a gamma- 
ray emitter with the High Energy Stereoscopic System (H.E.S.S.) array and the Fermi Large Area Telescope (Fenra'-LAT). 

Methods. Data obtained in 20.2 hours of dedicated H.E.S.S. observations and 38 months of Fenra'-LAT data, gathered by its usual all-sky scanning 
mode, have been analyzed to search for a gamma-ray signal. 

Results. No signal has been found in either data set. Upper limits on the gamma-ray flux are derived and are compared to models. These are the 
first limits on gamma-ray emission ever presented for galaxy clusters hosting cluster scale AGN outbursts. 

Conclusions. The non-detection of Hydra A in gamma-rays has important implications on the particle populations and physical conditions inside 
the bubbles in this system. For the case of bubbles mainly supported by hadronic cosmic rays, the most favorable scenario, that involves full mixing 
between cosmic rays and embedding medium, can be excluded. However, hadronic cosmic rays still remain a viable pressure support agent to 
sustain the bubbles against the thermal pressure of the ambient medium. The largest population of highly-energetic electrons which are relevant 
for inverse-Compton gamma-ray production is found in the youngest inner lobes of Hydra A. The limit on the inverse-Compton gamma-ray flux 
excludes a magnetic field below half of the equipartition value of 16 yuG in the inner lobes. 



1. Introduction 
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ated by these AGN on the embedding Intra-Cluster Medium 
(ICM) can be seen in se veral systems (for a review, see 
iMcNamara & Nulsenl 120071) . Typical signatures for an AGN - 
ICM interaction are surface brightness depressions in the dif- 
fuse thermal X-ray emission of the cluster which are caused 
by cavities in the ICM. These cavities appear to be filled with 
non-thermal electrons whi ch radiate in the radio band due to 
synch rotron emission (e.g. iBirzan et all 120041 : iDunn & Fabianl 
120061) . AGN-blown bubbles surrounded by thermal plasma offer 
the exciting possibility to constrain the energetics of these out- 
bursts. This can be done by estimating the work that is necessary 
to expand the bubbles against the thermal pressure of the embed- 
ding ICM (pV work in the following). The energetics involved 
in this AGN activit y can be enormous, i n som e cases even ex- 
ceeding 10 61 erg ( IMcNamara & Nulsenl l2007t) . The most pow- 
erful AGN outbursts kno wn to date are found in MS 0735+ 7421 
(IMcNamara et alll2005l), Hercu les A dNulsen et aUl2005al) and 
Hydra A dNulsen et aUl2005bl) . The AGN created bubbles in 
these systems have ages of about 10 8 years and exhibit sizes on 
the scale of the galaxy cluster itself. 

The nature of the main pressure support agent which fills 
the bubbles in the ICM is not known to date. Viable possibil- 
ities for the pressure support in such systems would be rela- 
tivistic particles such as hadronic cosmic ra ys or electrons (e.g . 
Dunn & Fabianl l2004t lOstrowski & SikoraL l200ll:lHinton et all 
2007), magnetic fields (e.g. lDunn & Fab ian. 2004) or hot plasma 



(e.g. lGittietaill2007h . The energy required to expand bubbles 
with volume V into a surrounding ICM with pressure p ranges 
from 2p V for magn etic fields to ApV for relativistic fluids such 
as cosmic rays (e.g. lWise et al l l2007h . 

One inevitable consequence of bubbles filled with non- 
thermal particles would be the production of gamma-ray emis- 
sion. For the case of hadronic cosmic rays, gamma rays are pro- 
duced by inelastic collisions between the high energy particle s 
and the thermal surrounding medium (e.g. iHinton et all 120071) . 
In case of electrons, gamma-rays are produced by up-scattering 
of cosmic microwave background (CMB) and infrared extra- 
galactic backgroun d light (EBL) photons by these electrons (e.g. 
lAbdo et ail l2010al) . Large-scale leptonic gamma-ray emission 
connected to AGN lobes has indeed been discovered with the 



Fermi satellite from the radio galaxy Centaurus A (Abdo et al 



and potentially from NGC 6251 (Ta keuchi et al 



Both galaxies are not hosted by a cluster. To date, 
no galaxy cluster has been firmly detected in gamma rays 
dPerkins et all 120061: lAharonian et all l2009allbl lAleksic et all 
l2010n Ackermann et al 1 120101) . However, the detection of an 
extended gamma-ray signal resulting from annihilation emis- 
sion from supersymmetric dark ma tter has be e n clai med for the 
Virgo, Fornax and Coma cluster (lHan et all |2012|) . Recently, 
NGC 1275 the central radio galaxy of the Perseus cluste r has 
been detected in VHE gamma rays (lAleksic et all 1201 2al) . and 
with this deep exposure stringent upper limits o n the emission 
of the Perseus cluster itself have been obtained (lAleksic et all 
l2012bl) . Galaxy clusters hosting cluster-scale AGN outbursts ap- 
pear to be promising targets for gamma-ray observations accord- 
ing to the extraordinary energetics inferred from the AGN - ICM 
interaction seen in these systems. To date, no gamma-ray obser- 
vations on galaxy clusters that host cluster-scale AGN outbursts 
have been presented. 

The Hydra A system (Abell 0780) at 
RA(J2000) = 9 h 18 m 05.7 s and Dec(J2000) = -12°05'44"at 
a redshift of 0.0538 is the closest kn own galaxy cluster wh ich 
hosts a cluster-scale AGN outburst dNulsen et all l2005bl) . It 
features several cavities with a total expansion work pV of 



4 x 10 60 erg done on the ICM. Thus the total energy required, 
depending on the equation of state of the main pressure agent, 
is (0.8 - 1.6) x 10 61 erg which were deposi t ed in the last few 
10 8 years in the surroundings dWise et a fl, 120071) . Hydra A 
also features low-freque ncy radio lobes ex tending to almost 4' 
from the cluster center dLane et all 20041). Sho cks in the ICM 
surround these radio lobes ( Nulsen et all 2005bl) with energetics 
of 9 x 10 60 erg, comparable to the expansion work done in the 
cavities against the thermal plasma. The central AGN outburst 
has also driven su bstantial gas dredge-up in the Hydra A system 
dGittiet alll20Tl . 

Chandra has furthermore revealed an extensive cavity sys- 
tem consisting of three generations of cavities with decreasing 
ages, which points towa rds a complex activity history of the sys- 
tem (Wise et a ll [2007b . Most relevant for gamma-ray produc- 
tion are the giant outer lobes that dominate the energetics in the 
Hydra A system and the inner lobes that are expected to contain 
the youngest population of particles. Both possibilities will be 
further discussed in Sec.|4] 

Due to its proximity and energetics, Hydra A is expected to 
feature the highest gamma-ray flux of all galaxy clusters har- 
boring cluster-scale AGN outbursts. For the case of hadron- 
dominated bubbles it was inferred that the flux might be close 
to the dete ction limit of the cur rent generation of gamma-ray in- 
struments dHinton et alll2007l) . 

In this paper, upper limits on the gamma-ray emission from 
the Hydra A system are reported. Limits obtained by the High 
Energy Stereoscopic System (H.E.S.S.) and the Fermi Large 
Area Telescope (ferm/-LAT) are presented in Sec. [2] and Sec. [3] 
respectively. These limits are used to obtain constraints on 
the energy of hadronic cosmic rays (Sec. 14.1. Il l and electrons 
(Sec. l4.1l2l and l4~2l i which may populate the AGN outburst re- 
gion in this galaxy cluster. 

Throughout this paper a ACDM cosmology with 
H = 70 km s _I Mpc" 1 , Q A = 0.7 and Q M = 0.3 is as- 
sumed, corresponding to a luminosity distance of = 240 
Mpc, an angular diameter dist ance of 216 Mpc and a linear 
scale of 1.05 kpc per arcsecond dWise et alll2007l) . 



2. H.E.S.S. data analysis 

Hydra A was observed in the VHE g amma-ray range 
(E > 100 GeV) with H.E.S.S. dHintonl \200A which is an ar- 
ray of four Imaging Atmospheric Cherenkov Telescopes located 
at the Khomas Highland in Namibia (23°16' 18" S 16°30'00" E, 
altitude 1800 m). Data were taken in March and April 2007 and 
from January to March 2010. In total 20.2 hours of good quality 
data (excluding data taken dur ing bad weather and da ta affected 
by hardware irregularities, see lAharonian et alll2006l) were col- 
lected. The data were obtained with a mean zenith angle of 15° 
which resulted in an energy threshold of 240 GeV. 

The data were analyzed with a boosted decision tree method 
(lOhm et alll2009l) . For H.E.S.S., Hydra A was treated as a point- 
like source. This is a reasonable assumption since the lobes 
of Hydra A extend over 4' in comparison to the 68% con- 
tai nment radius of t h e H.E .S.S. point spread function (PSF) of 
6' dAharonian et all 120061) . Using £ std cuts (£ denotes the 
boosted decision tree clas sifier and for the definition of std 
; Ohm et aLll2009in and reflected background model 



cuts seelL 

dBerge et al. . 20071) . a total number of counts Nq^ of 456 on 



the target and a number of background counts Noff of 7265 
(with source to background normalization a = 0.0614) were 



software version hap-1 l-02-pl07 
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measured. This results in an e xcess of 9 . 7 eve nts correspond- 
ing to a significance of OAcr dLi & MaL Il983h and hence no 
significant signal has been found (see Fig. [TJ. This result was 
confirmed with an independe nt calibration and analysis chain 
dde Naurois & Rollandl 12009*). Since no sig nal was detec t ed, up - 
per limits were derived using the method o f iRolke et al.l (120051) . 
At a confidence level of 95% an upper limit of F y (> 240 GeV) < 
7.9 x 10 -13 cm _2 s _1 for a power-law of the form dN/dE cc E~ r 
with an assumed photon index F = 2.5 is found. The gamma- 
ray index is chosen to approximately matc h the shape of the 
predicted spectrum (see lHinton et al.Ll2007l) . Upper limits were 
computed also for F = 2.0 and only weakly depend on the spec- 
tral index with a difference of less than 10%. In Fig. [2] upper 
limits for F = 2.5 that fits closest to the model predictions, are 
plotted and compared to model predictions. 



3. Fermi data analysis 

Hydra A has been observed by the Large Area Telescope (LAT), 
which is the primary instrument on the Fermi Gamma-ray Space 
Telescope (Fermi). It is a pair conversion telescope for high- 
energy gamma-rays with a wide field of view. It covers the en- 
ergy regime from 20 MeV to 300 GeV with an angular resolution 
of appro ximately 3.5° a t 100 M eV and narrowing to 0.14° at 10 
GeV (see lAtwoodelafl 120091) . In survey mode, the observatory 
is rocked north and south on alternate orbits so that every part of 
the sky is observed for ~30 minutes every 3 hours. 

Fermi-LAS' observations of Hydra A from MJD 54682.9 
to 55816.7, corresponding to a period of ~38 months from 
August 2008 to September 2011, are used in this paper. The 
data were retrieved from the public data archive and analyzed 
using the Fermi Science Tools v9r23 package. The standard 
event filt ering, recon s tructio n and classification were applied to 
the data dAbdo et all l2009h . The instrument response function, 
P7S0URCE_V6, is applied throughout the data analysis. 

Events with energies between 200 MeV and 200 GeV and 
within a circular region of 15° radius have been considered in the 
analysis. A binned maximum likelihood analysis was performed 
on the data using the gtlike t ool. A ll sources in the Fermi sec- 
ond year catalog dAbdo et al. Ll201lT) within 15° are included in 
the source model as well as the Galactic diffuse emission model, 
gal_2yearp7v6_v0 . fits, and the corresponding extragalactic 
isotropic diffuse emission model. The isotropic model used is 
iso_p7v6source . txt, which is valid for P7S0URCE_V6 instru- 
ment response functions. All point sources were modeled with 
parameters fixed to those reported in the Fermi second year cat- 
alog unless they were within 10° of Hydra A. Since Hydra A 
does not appear in the Fermi second year catalog, an additional 
point source was inserted at its position, assuming a power-law 
spectrum. Details on the likelihood analysis techniques and the 
models used can be found on the Fe rmi Science Support Center 
websitsQ (see also lAbdo etaill2009t) . 

The TS value, which square root is a measure of the signif- 
icance of a source, at the position of Hydra A is about 1 . This 
implies that no significant signal has been found in the data. 
Therefore, 95% flux upper limits are produced for a point-like 
source with an assumed spectral index of F (1.5, 2.0 and 2.5, 
that are given in Table [T]). Models are compared to upper limits 
obtained for the spectral index that is most compatible with the 
prediction. Flux upper limits are shown in Fig. [2] for F = 1.5 
(blue) and F = 2.0 (cyan) whereas the limit obtained for T = 2.5 
is shown in Fig [3] and |4] 



1.5 
2.0 
2.5 



F l ,i(>200MeV) [erg cnr^IrT 



1.4 x I0r u 

2.1 x lCr 13 

3.2 x irr 13 



Table 1. Fermi-LAT upper limits obtained for the AGN outburst 
in Hydra A. 
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Fig. 1. Significance map from the H.E.S.S. data . Overlaid are 
conto urs of the NVSS 1.4 GHz radio survey dCondon et all 
1 1998b which show the extension of the lobes in Hydra A. 
Additionally the 68% containment radius of the H.E.S.S. PSF 
is indicated by the dashed circle. 



4. Discussion 

Several scenarios of gamma-ray production in the Hydra A 
galaxy cluster are possible. The whole system consists of three 
generations of AGN outbursts that c reated bubbles of different 
energetics and consecutive ages (see lWise et aTll2007l) . The gi- 
ant outer lobes inflated by the oldest cycle of AGN activity con- 
tain most of the energy. For these outer lobes a hadronic scenario 
is investigated in Sec. 14.1.11 and a potential leptonic scenario is 
discussed in Sec. 14.1.21 In a hadronic scenario it is expected that 
gamma-ray emission is dominated by the giant outer lobes, since 
cooling of hadronic cosmic rays on timescales re levant for the 
evolu tion of the outer lobes of ~ 10 8 years (see iHinton et all 
l2007h is unimportant in the Hydra A system. The situation is 
different for electrons since they are effected by cooling on this 
rather long time scale. It is expected that the largest population 
of highly-energetic electrons which are relevant for IC gamma- 
ray production is found in the youngest inner lobes. Therefore, a 
leptonic scenario for the inner lobes is examined in Sec. 14.21 



http : //f ermi . gsfc . nasa . gov/ssc/ 



4.1. Constraints on the particle population in the giant outer 
lobes 

4.1.1. Hadronic population 

The gamma-ray luminosity in a hadronic scenario is given by the 
total energy in cosmic rays and the mean density of target ma- 
terial. In galaxy clusters which harbor cluster-scale AGN out- 
bursts, the total energy in cosmic rays can in principle be es- 
timated from the energetics of the AGN outbursts if this is as- 
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sumed to be the energetically dominant feedback agent on the 
1CM. Viable proxies for the energy in cosmic rays could be the 
energetics of the shock wave or the energy needed to sustain 
the X-ray cavities. In the general picture for AGN outbursts in 
galaxy clusters, the radio bubbles are dominated by cosmic rays, 
whereas the thermal ICM is distributed around these bubbles. 
The radio bubbles appear as surface brightness depressions in 
X-rays. This indicates a depletion of the hot ICM inside them. 
The three-dimensional structure of the bubbles is not known and, 
consequently, the actual density of X-ray emitting plasma in- 
side the cavities can only loosely be constrained. Limits on the 
density of thermal plasma inside the lobes can also be obtained 
with the depolarization effect of the radio emission of the lobes 
dGarrington & Conwavil200TI) . The actual gamma-ray luminos- 
ity of the system will depend on the level of mixing between 
hadronic cosmic rays and the thermal ICM. Processes which can 
lead to an effective mixing between cosmic rays and target mate- 
rial in AGN outbursts are diffusion of cosmic rays out of the bub- 
bles t o the regions with higher ICM density (e.g. lHinton et all 
l2007h . and entrainment of non-relativistic material in the outflow 
from the central engine (e.g. lPope et aLlUoiOl) . 

For Hydra A order of magnitude estimates give a gamma- 
ray luminosity of L y = £ pp /3t w ~ 10 43 erg s _1 , with E w ~ 
10 61 erg is the total energy in hadronic cosmic rays and t pp ~ 
6 x 10 9 years is the cooling time for proton-proton interac- 
tions for a mean density of target mate rial of 5 x 1CT 3 cm 3 , 
as obtained from X-ray measurements dNulsen et al.L l2005bl) . 
This results in a gamma-ray flux at Earth of F y = L y IAnd 2 L ~ 
10~ 12 erg cm~ 2 s _1 . This estimate shows that for an optimistic 
scenario where cosmic rays are well mixed with the embedding 
target material, this object is within reach of the current gener- 
ation of gamma-ray instruments. The predicted gamma-ray flux 
for a more elaborate model for Hydra A, assuming a hadronic 
emission mechanism and corrected for absorption by the EBL 
for a re dshift of 0.0538, is s hown for different scenarios (mod- 
els from lHinton et al.L 120071) in Fig. [2] For all cases the adopted 
mean density of thermal plasma outside the bubbles is 5 x 10~ 3 
crrT 3 . 

Cases (a) and (b) consider an energy in cosmic rays corre- 
sponding to the energetics of the blast wave surrounding the bub- 
bles of 9 x 10 60 erg. In case (a) the bubbles are filled with cold, 
unseen gas with the same density as the surroundings which 
could be entrained by the AGN outflow. In case (b) bubbles 
are completely evacuated from the thermal ICM and mixing 
between the thermal ICM, and the cosmic rays results solely 
from energy-dependent diffusion of cosmic rays to the outside 
medium of the bubbles. 

For the cases (c) and (d) it is adopted that the total energy in 
cosmic rays is IpV — 4 x 10 60 erg, which is necessary to prevent 
the cavities in the X-ray emitting gas to collapse. In case (c) the 
density of the ICM in the bubbles is half of the density outside 
the bubbles and for case (d) the same scenario comprising empty 
cavities as (b) is adopted. 

Predictions of these different model assumptions are com- 
pared to the upper limits on the gamma-ray emission of the 
Hydra A system obtained with Fermi-LAT and H.E.S.S. From 
Fig. |2] it is evident that these instruments are able to constrain 
the most favorable scenario for hadronic cosmic ray content and 
mixing of cosmic rays and thermal gas in the Hydra A galaxy 
cluster. This scenario would require a complete compound be- 
tween cosmic rays and ICM. This model (a) predicts a flux of 
E 2 dN/dE =s 4 x 10~ 13 erg crrT 2 s _1 in the range of about 1 GeV 
- 300 GeV. The presence of cavities in the ICM seems to ar- 



gue against complete mixing between these two components. 
However, it has to be noted that 10% of the entire ICM con- 
tained within a radius of 150 kpc from the cluster center has been 
dredged-up by the AGN outburst. This up-lifted coo ler gas par- 
tially follow the location of the giant outer bubbles dGitti et al.L 
|2012|) . Since these bubbles occupy 10% of the cluster volume 
within a radius of 150 kpc (Wi se et ail 120071) . significant en- 
trainment of cool gas in the outer bubbles can be expected. 
From the upper limits obtained with Fermi-LAT (see Tab.[TJ and 
H.E.S.S. above 240 GeV of 2.8 x 10~ 13 erg cirr 2 s" 1 (assuming 
F = 2.5) limits on the degree of mixing between cosmic rays 
and ICM can be derived. Fermi-LAT can constrain the degree 
of mixing to 0.5 and H.E.S.S. can limit the degree of mixing to 
less than 0.7 where means no mixing and 1 defines complete 
mixing between the two components. It is expected that parti- 
cles with higher energies diffuse faster into their surroundings 
and therefore also mix faster with the ambient medium. Thus 
both values for the limit on the degree of mixing at different en- 
ergies provide interesting constraints on particle transport in the 
Hydra A system. In general, hadronic cosmic rays as the ener- 
getically most important feedback agent in cluster-scale AGN 
outbursts can currently not be excluded. 
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Fig. 2. Upper limits on the gamma-ray flux are compared to the 
predicted spectral energy distribution for a hadronic scenario for 
the Hydra A system. Fermi limits are shown in blue for F = 1.5 
and cyan for F = 2.0 and H.E.S.S. limits are displayed in red for 
F = 2.5. Gamma-ray indices are chosen to approximately match 
the shape of the predicted spectrum. Upper limits for Fermi and 
H.E.S.S. for the entire probed energy range are shown with the 
assumed spectral index and for consistency the same represen- 
tation for the Fermi and H.E.S.S. limits are used. The contin- 
uation of the H.E.S.S. limit of 2.8 x 10~ 13 erg crrT 2 s _1 above 
the threshold of 240 GeV towards higher energies only reflects 
the adopted spectral index of T = 2.5 and therefore does not 
represent the general H.E.S.S. sensitivity at higher energies. For 
details of the different models see main text. Model (a) would 
predict an integral flux of F{> 240 GeV) * 1.5 x 10~ 12 crrrV 1 
above the H.E.S.S. threshold. 
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4.1.2. Electronic population 

Electrons, as opposed to hadrons, can lose their energy effi- 
ciently durin g the evolution time of the outburst in Hydra A of 
~ 10 8 years dHinton et all |2007|). For magnetic fields found in 
Hydra A of about 6 fj.G dTavlor & Perle"vl Tl993) only electrons 
with energies < 1 GeV are still present after synchrotron-cooling 
on such a time-scale. This picture is supported by a steepening 
of the radio index of synchrotron radiation from about a - 0.7 
close to the core to almost a 2 at the edge of the lobes in the 
radio band of 330 - 1415 MHz dLane et al.1 [2004 . For such a 
cooling-dominated electron population no gamma-ray inverse- 
Compton (IC) emission is expected in the Fermi-LAT energy 
range. 

However, Fermi-LAT has detected gamma-r ay emission 
from t he giant lobes of the Centaurus A system dAbdo et al.L 
[2010ah . Since this is the only system where gamma-ray emis- 
sion from an AGN outburst has been detected on spatial scales 
of the order of 100 kpc, it is the only example which can in prin- 
ciple be compared to the limits obtained for the Hydra A lobes. 
It has to be noted that Centaurus A and Hydra A are quite dif- 
ferent systems. In contrast to Centaurus A, Hydra A is located 
in a galaxy cluster environment with buoyantly rising bubbles. 
This fact together with different jet power, differing black hole 
mass and different accretion history in both systems may limit 
the applicability of such a comparison. 

The fermz'-LAT discovery of the Centaurus A lobes has been 
interpreted in the framework of a leptonic scenario where ener- 
getic electrons up-scatter CMB and EBL photons to gamma-ray 
energies. The limited radiative lifetimes of these energetic elec- 
trons may point towards in situ particle acceleration in the lobes. 
These processes may also be at work in the Hydra A system. 

To test a similar scenario for Hydra A as has been found for 
Centaurus A, the upper limits on gamma-ray emission obtained 
from Hydra A are compared to the flux measurements of the 
Centaurus A lobes scaled for the different distance and energet- 
ics of the Hydra A system. For this comparison for the combined 
emission of both Centaurus A lobes a flux above 100 MeV of 
1.86 x 10~ 7 ph cmy 2 s -1 and a spectral index of 2.55 is adopted 
dAbdo et al.L l2010al) . It has to be noted that in addition to the 
lobes, the nucleu s of Centaurus A als o emits gamma-rays with a 
comparable flux dAbdo et al l l20T0bl) but for the following dis- 
cussion only the emission from the lobes is considered. To con- 
strain the total energy in electrons in the lobes also the inten- 
sity of the photon field available for IC up-scattering has to be 
known. For the case of Centaurus A it was found that CMB and 
EBL are the dominating photon fields and a total energy of elec- 
trons of 1.5 x 10 58 erg was estimated dAbdo et al.Ll2010al) . The 
intensities of the CMB and EBL photon fields are equivalent for 
the Hydra A and the Centaurus A systems. With the assumption 
that the Hydra A lobes contain electrons with the same spectral 
characteristics as the Centaurus A lobes and corrected for the 
distance to Hydra A it is found that the total energy in electrons 
in the Hydra A lobes can be constrained to < 2 x 10 60 erg with 
an uncertainty dominated by the me asured gamma-ray fl ux of 
the Centaurus A lobes of about 30% dAbdo et al.L l2010al) . This 
is smaller than the 1 pV work of 4 x 10 60 erg (see Fig. [3}. This 
limit can be regarded as conservative since the central galaxy in 
Hydra A with log(Lv/[erg/s]) = 45.16 (apparent V-band mag- 
nitude my = 12.63, extinction Ay = 0.139^3 

is about 20 times 

more luminous than Centaurus A and can therefore provide an 
additional photon field for IC up-scattering (see Sec. l4.21 i. 



To summarize, for an electron population with the same 
spectral characteristics as it is found in the Centaurus A lobes, 
electrons can be excluded to be the main pressure support of the 
large scale bubbles in Hydra A. This result together with the fact 
that radio emission from the Centaurus A giant lob es extends to 
higher frequency than in in the Hydra A giant lobes dAbdo et all 
1201 Oat ILane et al.l |2004) points towards quite distinct properties 
of the electron population in these systems. 
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Fig. 3. Comparison between scaled gamma-ray flux of the 
Centaurus A lobes and the upper limits for Hydra A. Fermi and 
H.E.S.S. limits are shown for F = 2.5. Gamma-ray indices are 
again chosen to approximately match the shape of the predicted 
spectrum. Here it is assumed that the Hydra A lobes contain elec- 
trons with the same spectral characteristics as the Centaurus A 
lobes and that CMB and EBL are the dominating photon fields 
which are up-scattered to gamma-ray energies. The gamma-ray 
flux of the Centaurus A lobes is scaled to a distance of 240 Mpc 
and to a total energy in electrons of 10 60 erg (solid line) and 
1 pV = 4 x 10 60 erg (dotted line), respectively, to account for the 
different distance and energetics of the Hydra A system. 



4.2. Leptonic scenario for the inner lobes 

The situation for leptonic gamma-ray emission may be more 
promising in the youngest, innermost lobes. Indeed the Hydra A 
system is detected in the radio band up to f requencies of 90 GHz 
dCotton et al.L I2009L IWright et all 120091) . These observations 
show the presence of highly-energetic electrons which are rel- 
evant for IC gamma-ray production. Thus, here the upper limits 
on the gamma-ray emission for Hydra A are compared to the 
expected gamma-ray flux from the inner radio lobes. 

For estimating the expected gamma-ray IC luminosity in a 
leptonic scenario, first the population of electrons in the lobes 
is explored according to their radio synchrotron emission. The 
particle populations are evaluated separately for the inner and 
o uter lobes of Hydr a A. The inner lobes (denoted as 'A' and 'B' 
in I Wise et al.L[2Q07h are assumed to be prolate spheroids located 
at the distance of r; n = 25 kpc from the center, with semi-minor 
axis a m = 20" and semi -major axis b \ „ = 35 " (as inferred from 
1.4 GHz radio maps in Birzanetal., 2008). The total volume 
of the inner lobes is therefore Vj n = 2 x (4/3)7ra?J?j n . For the 
outer lobes (denoted as 'E' and "F' in lWise et all 20071) r^ , t = 
225 kpc, fl out = 90" and b out = 120" dBirzan et all 120081) are 
used. 
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The rad i o flux e s 
IBirzan et al.l (I2008h 



for the inner lob es are taken from 
ICotton et all (120091) . and IWright et al.l 



d2009l) and are shown as red open circles, red filled circles, and 
red stars in Fig. [4] r espective l y. The radio flux for the outer lobes 
is adopted from lBirzan et al.ld2008l black square i n Fig .|4]i. Note 



that the WMAP fluxes as given in IWright et"aH d2009l) match 
well the radio continuum of the inner lobes, with the exception 
of the 61 GHz flux. Therefore it is assumed that these fluxes rep- 
resent indeed the high-energy tail of the synchrotron emission of 
the inner lobes. For the 6 1 GHz flux as well as the 90 GHz flux 
dCotton et al.l 120091) it is anticipated that they are dominated by 
the fiat-spectrum synchrotron emission of the jets and the nu- 
cleus of the radio g alaxy in Hydra A in stead of the inner lobes. 

As discussed in lBirzan et al.l(l2008l) the synchrotron continua 
of both the inner and the outer lobes of Hydra A can be well rep- 
resented by broken power laws with the low- and high-frequency 
radio spectral indices ai ow =; 0.5, and ahi g h - 1-5, respectively. 
Therefore the lobes' electron energy distribution is assumed to 
be n e {y) oc y~ 2 for 1 < y < y bl , and n e {y) oc y~ 4 x exp[-y/y max ] 
for y > ybi ■■ In the case of the inner lobes, the break and the max- 
imum electron Lorentz factors are assumed to correspond to the 
synchrotron frequencies 10 GHz and 100 GHz, respectively. In 
the case of the outer lobes the analogous synchrotron break and 
maximum frequencies are taken as 0.1 GHz and 1 GHz, respec- 
tively. 

For the evaluation of the IC-fluxes the following target pho- 
ton fields are taken into account: CMB, EBL, IR emission of the 
nuclear dust, and the starlight emission of the elli ptical host. The 
EBL i s approximated by the spectrum given in iRaue & Mazinl 
d2008l) . For the circumnuclear dust emission a modified black 
body spectrum with the dust temperature 60 K and the total IR 
luminosity integrated over the frequency range 10 11 - 10° Hz 
equal to Ldust — 7 X 10 43 erg s _1 is assumed. This model can ac- 
count well f or the SCUBA and MIPS data for the radio galaxy 



inHvdra A (IShi et all 120051 IZemcov et"aIL 120071: iDicken et all 
120081) . The IR energy density at the position of the inner and 
outer lobes is therefore taken as t/dust = ^dust/47ir 2 c. In the case 
of the starlight emissi on of the host galaxy , the template spec- 
trum as discussed in IStawarz et al.l (120061) . normalized to the 
total V-band luminosity Ly = 1.45 x 10 45 ergs~', is adopted. 
The energy density of the starlight at the position of the lobes 
is then evaluated as U stm - = L star /4^r 2 c. For the inner lobes the 
dominating photon fields are CMB in the microwave regime, the 
dust emission in the infrared and the star light in the optical, 
respectively. EBL has been found to be unimportant for IC up- 
scattering in the inner lobes. 

With all the model parameters and assumptions as specified 
above, the synchrotron and inverse-Compton fluxes of the in- 
ner and outer lobes of Hydra A are evaluated for the remaining 
two free parameters: the lobes' magnetic field intensity B, and 
the equipartition ratio tj. The latter parameter is defined as r\ = 
U c /Ub where U e = J dyym e c 2 n e (y), and Ub = B 2 /8n. These 
values also define the total energy in ultrarelativistic electrons 
and the magnetic field £ e +B = V x (U B + U e ) = V x t/ B x (1 + 77) 
stored in the bubbles with volume V. The calculations are done 
for different sets of the values of B and 77 to match the observed 
radio fluxes in all the cases. The results are presented in Fig. [4] 
There the red curves correspond to the emission of the inner 
lobes, and the black curves represent the outer lobes. The eval- 
uated synchrotron fluxes (see Fig. [4]i match well the collected 
data-set and ar e in genera l agree ment with the spectral an alysis 
carried out by lLane et all d2004l) and iBirzan et all d2008l) . The 
solid curves illustrate the case with the exact electron-magnetic 
field energy equipartition, namely B{ n = I6/1G, rji n = 1, B ollt = 
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Fig. 4. Synchrotron and expected corresponding IC emission for 
the inner and outer lobes of H ydra A. The experimental data for 
the inn er lo bes are taken from Birzan et"ail d2008l) . ICotton et al.l 
(120091) . and lWrightetail d2009i and are shown as red open cir- 
cles, red filled circles, and red stars, respectively. The radio flux 
for the outer lobes (indicated by black squares) is adopted from 
IBirzan et al.l (|2008). Additionally the upper limit for the power- 
law emission at 1 keV for the outer (Northern) lobe (black arrow; 
lHardcastle & Crostonl 120101) is included. The red curves corre- 
spond to the emission of the inner lobes, and the black curves 
represent the outer lobes. The solid curves illustrate model (a), 
the dashed curve model (b) and the dotted curve is for model (c) 
from Tab. [2] Fermi and H.E.S.S. limits are shown for F = 2.5. 



6.3 //G, and j] out = 1 (model a). Note that the standard minimum- 
energy calculations return typically the equipartition magnetic 
field 15 - 3 /uG for the inner l obes and 3 - 6/uG for the outer 
lobes (e.g., iTavlor et all 1 19901: IBirzan et all 120081) which is in 
agreement with our modeling. Parameters of further models can 
be found in Tab. [2] Model (b) is for a case with the magnetic 
field twice lower than the equipartition value, and in model (c) 
the magnetic field is three times below the equipartition value. 



model 




7/in 






^/out 


£e+B,out 




UiG] 




[erg] 






[erg] 


a 


16 


1 


4.1 x 10 58 


6.3 


1 


4.3 x 10 59 


b 


8 


12 


1.2 x 10 59 


3 


12 


1.2 x 10 60 


c 


5 


65 


2.6 x 10 59 


2 


65 


2.9 x 10 60 



Table 2. Properties of the models to calculate the IC emission 
from the inner and outer lobes. B gives the magnetic field inside 
the lobes, 77 is the equipartition ratio and E e+ B is the total energy 
in electrons and magentic field for the inner (subscript in) and 
outer (subscript out) lobes, respectively. 



These calculations are now compared to the upper limits ob- 
tained in the gamma-ray range. Additionally also the upper limit 
for the power-law emissi on at 1 keV for the oute r (Nort hern) lobe 
(black arrow in Fig. H lHardcastle & Crostonl l2010h obtained 
with XMM is included. The inverse-Compton emission of the 
lobes in Hydra A is expected to be negligible at TeV photon en- 
ergies. It may however be pronounced within the lower energy 
range from keV-to-GeV photons. The XMM and ferm/-LAT 
upper limits for the Hydra A system seem to already exclude a 
magnetic field in the lobes below half of the equipartition value. 
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5. Conclusion 

In this paper the nearby galaxy cluster Hydra A that hosts a 
cluster-scale AGN outburst is investigated for being a gamma- 
ray emitter. Galaxy clusters hosting a cluster-scale AGN outburst 
are potentially detectable gamma-ray sources due to the enor- 
mous energetics inferred from the observed AGN - ICM inter- 
actions. However, only upper limits could be obtained from 20.2 
hours of H.E.S.S. and 38 month of Fermi-LAT observations. The 
non-detection of Hydra A in gamma-rays has important implica- 
tions on the particle populations and physical conditions inside 
the bubbles in this system. These upper limits constrain the total 
energy contained in relativistic particles such as hadronic cosmic 
rays and electrons, which can be compared to the energy which 
is necessary to prevent the observed cavities in the ICM from 
collapsing. 

Constraints on the particle population in the Hydra A 
galaxy cluster can also be compared to the limits on such 
a non-thermal comp onent inferred for the Perseus cluster 
( lAleksic et al.l l2.Q12bh . The Perseus cluster is an interesting can- 
didate for this comparison since it also shows signatures of AGN 
- IC M interaction s in for m of radio lobes and cavit ies in the ICM 
e.g. Fabian et al., For the Perseus cluster Aleksic et al. 

2012bl constrained the average fraction of energy in hadronic 
cosmic rays to thermal energy Eqr/E^ to < 1 - 2% depend- 
ing on the exact assumptions. For Hydra A the hadronic cos- 
mic ray content in the central 200 kpc can be limited to about 
5 x 10 60 erg, assuming complete mixing between cosmic rays and 
ICM (see Fig.|2|. When adopting a total gas mass of 5 x 1O 12 M 
and a temperature of 4 keV for the central 200 kpc of Hydra A 
(iDavis et all 1200 lh then Ec R /E th < 13% is found. This is signif- 
icantly less constraining than for the case of the Perseus cluster. 
Hydra A, however, is the prime candidate to explore the particle 
content of giant AGN-blown lobes in galaxy clusters. This fol- 
lows from the fact that the AGN outburst in the Perseus cluster 
is about an o r der of magnitude less energetic (pV » 3 x 10 59 erg, 
iFabian et al.L 1201 lh than the AGN feedback in Hydra A. The 
smaller energetics is not readily compensated by the shorter dis- 
tance to the Perseus cluster (c/l = 75 Mpc). 

For Hydra A for the case of bubbles mainly supported by 
hadronic cosmic rays these upper limits can exclude the most 
favorable model, that requires full mixing between relativistic 
particles and embedding thermal medium. It is found that Fermi- 
LAT can constrain the degree of mixing to 50% and H.E.S.S. can 
limit the degree of mixing to less than 70%. However, hadronic 
cosmic rays still remain a viable pressure support for the bub- 
bles. 

In contrast to hadrons, electrons cool quite fast above 
GeV energies in the environment of the Hydra A system. 
Consequently, a passively evolving population of electrons in the 
oldest outer lobes cannot be detected with the presented observa- 
tions. However, for the youngest, inner radio lobes, the limit on 
the IC flux seems to exclude a magnetic field below about 8 fiG, 
that is half of the equipartition value. For the large outer lobes, 
a population of electrons rejuvenated by in situ particle acceler- 
ation comparable to the one detected in the Centaurus A system 
can be excluded as the main pressure support of the bubbles. 
Upper limits in the VHE gamma-ray range are not constraining 
for leptonic scenarios with respect to limits obtained in the GeV 
range. 

The main feedback agent which drives the evolution of the 
cavities in the ICM in the Hydra A galaxy cluster still remains 
unidentified. The upcoming Cherenkov Telescope Array (CTA, 
lActis et ail 1201 lh with its increased sensitivity will be crucial 



to test especially the presence of hadronic cosmic rays in the 
Hydra A galaxy cluster. 
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